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p ■ Abstract 

^ ; Anticipating the important role of tree level FCNC processes in the indirect search for new 
physics at distance scales as short as 1CT 19 — 1CT 21 m, we present complete NLO QCD 
corrections to tree level AF = 2 processes mediated by heavy colourless gauge bosons 
{Sj [ and scalars. Such contributions can be present at the fundamental level when the GIM 
mechanism is absent as in numerous Z' models, gauged flavour models with new very 
heavy neutral gauge bosons and Left-Right symmetric models with heavy neutral scalars. 
They can also be generated at one loop in models having GIM at the fundamental level 
and Minimal Flavour Violation of which Two-Higgs Doublet models with and without su- 
persymmetry are the best known examples. In models containing vectorial heavy fermions 
that mix with the standard chiral quarks and models in which Z° and SM neutral Higgs 
H° mix with new heavy gauge bosons and scalars in the process of electroweak symmetry 
breaking also tree-level Z° and SM neutral Higgs H° contributions to AF = 2 processes 
are possible. In all these extensions new local operators absent in the SM are generated 
having Wilson coefficients that are generally much stronger affected by renormalization 
group QCD effects than it is the case of the SM operators. The new aspect of our work 
is the calculation of 0(a s ) corrections to matching conditions for the Wilson coefficients 
of the contributing operators in the NDR-MS scheme that can be used in all models 
listed above. This allows to reduce certain unphysical scale and renormalization scheme 
dependences in the existing NLO calculations. We show explicitly how our results can be 
combined with the analytic formulae for the so-called P° QCD factors that include both 
hadronic matrix elements of contributing operators and renormalization group evolution 
from high energy to low energy scales. For the masses of heavy gauge bosons and scalars 
(9(1) TeV the remaining unphysical scale dependences for the mixing amplitudes M 12 are 
reduced typically from 10 — 25%, depending on the operator considered, down to 1 — 2%. 
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1 Introduction 

The absence of tree level flavour changing neutral currents (FCNC) within the Standard 
Model (SM), known under the name of the GIM mechanism [I], played a very important 
role in the construction of this model and undoubtedly contributed to its successes in an 
important manner. Not only are tree level FCNC processes mediated by Z boson and 
the Higgs absent in this model, but also the breakdown of the GIM mechanism at the 
one-loop level, governed by the hierarchical structure of quark masses and of the CKM 
matrix, appears to be an adequate description of existing data on FCNC processes within 
the experimental and theoretical uncertainties. 

Beyond the SM the GIM mechanism ceases to be a general property and there exist a 
number of popular models in which FCNC processes take place already at the tree level. 
The best known are various versions of the so-called Z' models [2] in which new neutral 
heavy weak boson (Z') mediate FCNC processes already at tree level. Similarly the heavy 
gauge bosons in the gauged flavour models [3] imply FCNCs at the tree-level. This is also 
the case in models based on left-right symmetry, where tree level heavy neutral Higgs 
exchanges contribute to AF = 2 amplitudes. 

Effective tree-level contributions to AF = 2 observables can also be generated at one 
loop in models having GIM at the fundamental level and Minimal Flavour Violation of 
which Two-Higgs Doublet Models with and without supersymmetry are the best known 
examples. In models with heavy vectorial fermions that mix with the standard chiral 
quarks and models in which Z° and SM neutral Higgs H° mix with new heavy gauge 
bosons and scalars in the process of electroweak symmetry breaking also tree-level Z° 
and SM neutral Higgs H° contributions to AF = 2 processes are possible. In all these 
extensions new local operators absent in the SM are generated having Wilson coefficients 
that are generally much stronger affected by renormalization group QCD effects than it 
is the case of the SM operators. 

Extensive model independent analyses of various authors of AF = 2 processes, in partic- 
ular in jU \5\ E], demonstrate very clearly that in the presence of 0(1) FCNC couplings 
of heavy gauge bosons and scalars, their masses must be above 100 TeV, correspond- 
ing to distance scales as short as 10~ 21 m, in order to satisfy the present experimental 
bounds. With couplings significantly suppressed, these masses can be lowered to the 
1 TeV — 10 TeV range and even lower scales, which are in the LHC reach. 

While until now no definite signs of new physics have been observed at the LHC, we expect 
that in the coming years new phenomena, new particles and forces, will be discovered and 
their nature tested both in high energy processes governed by ATLAS and CMS and 
low-energy high precision experiments with prominent role played by LHCb, Belle II, 
K + — > -k + vv experiment at CERN and later by the Super-B factory in Rome and the 
X-project at Fermilab. 

It is conceivable that these experiments will answer some of the present questions, si- 
multaneously opening new ones that will require to search for new physics beyond the 
reach of the LHC. While loop diagrams, like penguin diagrams of various sorts and box 
diagrams dominated the physics of flavour changing neutral current (FCNC) processes in 
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the last thirty years both within the SM and several of its extensions, we should hope 
that in the case of new particles with masses above 10 TeV this role will be taken over 
by tree-level diagrams. The reason is simple. Internal particles with such large masses, 
if hidden in loop diagrams, will quite generally imply very small effects that will be very 
difficult to measure. On the other hand tree diagrams could still provide a large window 
to these very short distance scales. 

Anticipating this future role of tree level diagrams we make another look at the NLO 
QCD corrections to AF = 2 processes like K° — K° and B® s — B® s mixings mediated by 
tree level heavy neutral gauge bosons and scalars. New contributions of this type imply 
the presence of new four-fermion operators in addition to the SM (V — A) x (V — A) 
operator. They have been classified in [Tl El El ID] and we will list them below in the basis 

of pig. 

Concerning QCD corrections, what is known are the renormalization group evolution 
matrices at the NLO level and the values of the hadronic matrix elements calculated 
using lattice methods. This information allows to study QCD effects in AF = 2 processes 
in a meaningful manner because only at the NLO level the Wilson coefficients can be 
properly combined with the hadronic matrix elements calculated by lattice methods at 
low energy scales. 

Now as pointed out in pT] instead of evaluating the hadronic matrix elements at the low 
energy scale we can choose to evaluate them at the high scale /Xj n at which heavy particles 
are integrated out. Thus the amplitude for a given M — M mixing (M = K, B^, B s ) is 
given simply by 

A(M^M) = (M\H a f= 2 \M) =Kj2cn^)(M\Q^^ m )\M), (1) 

i.a 

with k specified below. Here the sum runs over all the contributed operators which will 
be listed in Section [2j The matrix elements for B^ — B% mixing are for instance given 
then as follows pH [12] 

(B^Q^M) = ± mBd F%P?{B d ) = P^in) , (2) 

where the coefficients P^(B d ), in which the scale /z in has been suppressed, collect com- 
pactly all RG effects from scales below n in as well as hadronic matrix elements obtained 
by lattice methods at low energy scales. Analytic formulae for these coefficients are given 
in [11] while the recent application of this method can be found in [131 EH E3 ESJ EZ] • As 
the Wilson coefficients Ci(fi in ) depend directly on the loop functions, tree diagram results 
and fundamental parameters of a given theory, this formulation is very transparent and 
interesting short distance NP effects are not hidden by complicated QCD effects. 

In this approach the hadronic matrix elements in Eq. (JTJ) are usually calculated in a 
particular renormalization scheme at the matching scale /i; n . This scale, while being of 
the order of the masses of heavy particles that have been integrated out, does not have 
to be equal to these masses. As the amplitude on the l.h.s of Eq. (TI}) cannot depend 
on the choice of the renormalization scheme and on the precise value of the scale /ii n , 
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these unphysical dependences have to be cancelled by the ones present also in the Wilson 
coefficients Cf(fi m ). To this end these coefficients have to be known at the NLO level 
which requires the calculation of 0(a s ) corrections to penguin diagrams, box diagrams 
and in particular tree diagrams in the full theory and matching this calculation to the 
corresponding effective theory. 

Now in most applications to date the coefficients C^{^- in ) in the extensions of the SM have 
been calculated at the leading order, leaving some left-over unphysical scheme and scale 
dependences in the resulting physical amplitudes. While presently these dependences 
are significantly smaller than the present uncertainties in the evaluation of the hadronic 
matrix elements, the situation could change in this decade. In the SM these coefficients 
are known at the NLO and in a few processes at the NNLO level. An up-to-date review 
can be found in |18j . 

The goal of our paper is the evaluation of Cf(n m ) for the cases of tree-level colourless 
neutral gauge boson and neutral scalar exchanges at the NLO level. This amounts to 
the calculation of 0(a s ) corrections to the tree diagrams in question. We will also show 
explicitly how our results should be combined with the coefficients P" so that our final 
results will be mixing amplitudes at the NLO level that are general enough to be used for 
any model in which tree-level contributions to AF = 2 processes from colourless neutral 
gauge boson and scalar exchanges are present. The analysis of coloured gauge bosons 
and scalars is in progress but the calculations in this case are more involved and we will 
present them in a separate publication [T9] . 

Our paper is organized as follows: In Section [2] we recall the general structure of the 
effective Hamiltonian for AF = 2 processes and we give the full list of four-fermion 
operators that contribute to these transitions. In Section |3j the most important section 
of our paper, we describe the calculation of one-loop QCD corrections to the coefficients 
in the NDR-MS scheme and collect our results. The separate results for the amplitudes 
in the full theory and the effective theory given in the Appendices C and D should enable 
interested readers to check our calculations. In Section H] we demonstrate analytically that 
the new contributions remove unphysical scale dependences in the formulae present in the 
literature. In Section we combine our results with the P" QCD factors of [TJ] obtaining 
in this manner the complete NLO results for the AF = 2 amplitudes governed by tree- 
level exchanges of colourless gauge bosons and scalars. This section can be considered as a 
compendium of the master formulae for tree-level contributions to the mixing amplitudes 
at the NLO level that are valid in any extension of the SM in which such contributions are 
present. In Section[6]we demonstrate numerically that the resulting amplitudes practically 
do not depend on the detail choice of the matching scale. We conclude with a brief 
summary in Section [7J 
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2 Theoretical Framework 
2.1 Preliminaries 

While in the SM only one operator contributes to each AF = 2 transition in the K 
and B systems, in the tree level FCNC transition considered here there are eight such 
operators of dimension six. Consequently the renormalization group (RG) QCD analysis 
becomes more involved and due to the presence of right-handed and scalar currents and the 
resulting structure of the new operators QCD corrections play a much more important role 
in new physics contributions than in the SM contributions. Therefore also the unphysical 
renormalization scheme and scale dependences are much more pronounced when not all 
NLO QCD corrections are taken into account. 

In what follows, after listing all contributing operators we will summarize the effective 
Hamiltonian for AF = 2 transitions. 



2.2 Local Operators 

The contributing four-fermion operators can be split into five separate sectors, according 
to the chirality of the quark fields they contain. For definiteness, we shall consider opera- 
tors responsible for the K°-K° mixing. The operators belonging to the first three sectors 
(VLL, LR and SLL) read [TO] : 



Q\\K) 



{s a P L d a )(s?PndP), 



(3) 



(s a P L d a )(sPp L d?), 

(s^PLd^a^P^) 

where = 5(7^, 7^] and Pl,r = §(l~F7s)- The operators belonging to the two remaining 
sectors (VRR and SRR) are obtained from QY LL and Qf LL by interchanging Pl and Pr. 
In the SM only the operator QY LL (K) is present. The operators relevant for B q (q = d, s) 
are obtained by replacing in (JHJ) K by B q , s by b and d by q. 



Qi(K) = 
Ql hh (K) = 



2.3 Effective Hamiltonian 

The effective Hamiltonian for AF = 2 transitions can be written in a general form as 
follows 

n A e t 2 = kJ2c?(»)Q? , (4) 

i,a 
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where Q\ are the operators given in Eq. (jBJ) and Cf (/i) their Wilson coefficients evaluated 
at a scale // at which the hadronic matrix elements are evaluated. The overall factor 
k depends on the contributing particles and will be chosen such that for non-vanishing 
Wilson coefficients Cj a (/i in ) = 1 in the LO. The scale \i can be the low energy scale 
fiL at which actual lattice calculations are performed or any other scale, in particular the 
matching scale ii- m as in Eq. (CQ). In this case the matrix elements are obtained by evolving 
by means of RG equations the lattice results from /xz, to /i; n . The result of this evolution 
is given in Eq. ([2j). The general NLO formulae for the coefficients as functions of the 
QCD coupling constant and the hadronic E>i parameters calculated by lattice methods 
are presented in [TT] . 

As already emphasized in the Introduction, this formulation is very powerful as it applies 
to any extension of the SM. What distinguishes various NP scenarios are 

• the contributing operators 

• their Wilson coefficients C 4 a (/ii n ) which depend directly on the fundamental param- 
eters of a given theory. 

The results of tree-level and loop calculations performed at the matching scale fi- m at 
which the heavy particles are integrated out depend explicitly on the these fundamental 
parameters which allows to see very transparently the short distance NP effects that are 
not hidden by complicated QCD effects which necessarily take place between high energy 
and low energy scales. 

2.4 The Operator Structure from Tree Level Exchanges 

In the present paper we will consider FCNC AF = 2 amplitudes generated through 
tree-level very heavy gauge boson and scalar exchanges independently whether flavour 
violating neutral couplings in question have been generated at the fundamental level or 
through loop corrections. The only assumption that we will make in the present paper is 
that exchanged neutral gauge bosons and scalars are colourless as in many NP scenarios 
listed above. 

It is instructive to compare the operator structures in the effective Hamiltonian for AF = 2 
transitions at scales 0(/j,- m ) resulting from tree- level exchanges that differ when gauge 
bosons and scalars with colour or without colour are exchanged. We have: 

• A tree level exchange of a colourless gauge boson with LH and RH couplings gener- 
ates at /i in the operators QY Lh , Qi RR and Q\ R . This is an example of Z' models and 
gauge flavour models [3J. Also tree- level Z° can be generated in some extensions 
of the SM, in particular when new heavy neutral gauge bosons mix with Z° and 
heavy vectorial fermions mix with chiral quarks of the SM. When QCD corrections 
at the taken into account also the operator Q\ K is generated but its Wilson 
coefficient is suppressed by a s (fj,- m ) relative to other operators as we will see in the 
next section. 
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• A tree level exchange of a gauge boson carrying colour generates the operators QY LL , 
QY rr , Q\ r and Q^ R even without including QCD corrections. An example is the 
tree-level exchange of the KK-gluon in the RS models. 

• A tree level exchange of a colourless Higgs scalar generates the operators Qf LL , Qf RR 
and Q\ R . When QCD corrections at the taken into account also the operators 
Q2 LL , Q| RR and Q\ R are generated but their Wilson coefficient are suppressed by 
a s (//j n ) relative to other operators as we will see in the next section. Such heavy 
scalars are present in supersymmetric models and in left-right symmetric models. 
Tree level exchanges of the SM H° could also be generated in certain models. 

• A tree level exchange of a Higgs scalar carrying colour generates the operators , 
Qf RR and Q\ r even without the inclusion of QCD corrections. 

As already stated before we concentrate here on the colourless gauge bosons and scalars. 
The case of coloured gauge bosons and scalars will be discussed elsewhere [T9] . 



3 Matching Conditions 
3.1 Preliminaries 

The calculations of 0(a s ) QCD corrections to Wilson coefficients are by now standard and 
have been described in several papers. In particular in Section 5.4.2 of [20] all necessary 
steps have been presented in detail in the case of charged currents within the SM, while 
[21] presents the calculation of 0(a s ) corrections to C^ LL within the SM. The novel feature 
of the calculations present below when compared with these two papers is the appearance 
of new operators but the procedure is the same: 

Step 1 

We first calculated the amplitudes in the full theory. They are given in the case of a gauge 
boson exchange and a scalar exchange in Figs. [1] and [2, respectively. In the presence of 
massless gluons one encounters infrared divergences. We have regulated these divergences 
by a common external momentum p with — p 2 > for all external massless fields as done 
in [20]. Equally well they could be regulated by setting all external momenta to zero but 
giving the external quarks non- vanishing masses as done in [21] . As the Wilson coefficients 
cannot depend on the employed infrared regulator, the same result should be obtained 
in both cases. In fact in the case of gauge boson exchanges we have performed also 
calculations with the mass regulator obtaining the same results for the Wilson coefficients. 

The ultraviolet divergences present in the vertex diagrams in Figs. Q] and [2] have been 
regulated using dimensional regularization with anti-commuting 75 in 4 — 2e dimensions. 

Step 2 

We have calculated the matrix elements of contributing operators by evaluating the di- 
agrams in Fig. [3] making the same assumptions about the external fields as in the first 
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step. In contrast to step 1 one has to renormalize the operators. This we do in the MS 
renormalization scheme with anti-commuting 75, which corresponds to the NDR scheme 
of [22] used also in [Hi] and pi]. 

Step 3 

We finally inserted the results of the two steps above into the formula like the one in 
Eq. ([1]) and comparing the coefficients of operators appearing on the l.h.s (full theory) 
and r.h.s (effective theory) we found the coefficients C?(yUi n ). As these coefficients cannot 
depend on the infrared behaviour of the theory, the dependences on p 2 found in the first 
two steps have to cancel each other in the evaluation of Cf (fj,- m ). Indeed we verified this 
explicitly. The interested reader can do this as well by inspecting our intermediate results 
that we present in Appendices C and D. The appearance of the renormalization scale ^ n 
can be traced to the use of dimensional regularization and the renormalization in the MS 
scheme. 

Very often in analyses in which NP contributions are governed by box diagrams the overall 
factor in front of the sum in Eq. ([1]) is chosen as in the SM. However, in our analysis it will 
be more convenient to use in each case the normalization in which the Wilson coefficient 
of the leading operator evaluated at the matching scale is equal to unity in the absence of 
QCD corrections. In this manner the applications of our formulae in various models will 
be facilitated. In what follows we will first present the general structure of the effective 
Hamiltonian in each case. Subsequently we will list our results for the Wilson coefficients 
including 0(a s ) corrections. 




Figure 1: Tree level diagram and one-loop QCD corrections to AS = 2 transition mediated 
by a gauge boson in the full theory. The mirror diagrams are not shown. 
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Figure 2: Tree level diagram and one-loop QCD corrections to AF = 2 transition mediated 
by a scalar particle in the full theory. The mirror diagrams are not shown. 




Figure 3: Leading order and one-loop diagrams in the effective theory. The mirror dia- 
grams are not shown. 

3.2 Results 

Defining the general couplings through the Feynman rules in Fig. H] we find the following 
results. 



3 Matching Conditions 



9 



i^Saf, [A l l(A)P L + A%(A)P R ] 

3p 




H iSap[A^{H)P L + A%{H)P R ] 

Figure 4: Feynman rules for colourless neutral gauge boson A with mass Ma, and neutral 
colourless scalar particle H with mass M H , where i, j denote different quark flavours and 
a, (3 the colours. 

3.2.1 Colourless gauge boson 



^r 2 =^g^^ VLL (/^r + ^g^r K (/^r H 

(5) 



Ml 



The operator Q\ K is only generated at one-loop but not present at tree level. Consequently 
its Wilson coefficient is 0(a s ). We find for an arbitrary number of colours N 



C - M = ™ = l + ^(l-i) (-3.o g f + H 

a s ( , Ml 11 

t t> a s 3 / , M\ 1 

^^ 1 + 4^iv(- l0g ^-6 



(6) 



a„ / . M\ _ 1 
4tt V "° A* 2 6 



(7) 



1 + -^ -log 



^) = ^(-61og^-l). (8) 
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3.2.2 Colourless scalar 



(Af(#)) 2 



2M 2 H 



[C SRR (iu)Ql RR + Cl RR (fi)Q s 2 RR ] (9) 



The operators Q2 LL//SRR and Q\ R are generated by QCD corrections. We find for an 
arbitrary number of colours N 



C 1 SL ^) = C 1 SR V) = l + ^(-31og^ + ^) , (10) 



a s 2-N 
47 AN 



/ Ml 3\ 



a s ( 1 , M 2 1 



(11) 



c 1 LR M = -|S- < 12 > 

The formulae presented in this subsection are the main results of our paper. 



4 Renormalization Scale Dependence 

One of the main virtues of our calculation of 0{a s ) corrections to Wilson coefficients 
at the high energy matching scale ji in is the cancellation of the jj,- in dependence of the 
renormalization group evolution matrix by the /i; n dependence of the Wilson coefficients in 
question. This cancellation requires particular values of the coefficients of the log(M 2 //if n ) 
in Cj(/ii n ) where M stands for the mass of a heavy gauge boson or heavy scalar involved. 
As this cancellation constitutes an important test of our results it is useful to derive a 
general condition on the coefficients of log(M 2 //i 2 n ) in Cj(// in ). 

To this end let us look as an example at the evolution matrix U (fib, fJ>m) defined by 

d(n b ) = U(n b , At in )C(// in ), (14) 
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Here C (/i) is a column vector of Wilson coefficients. Expanding then this matrix around 
the two fixed scales m& and M keeping only the logarithmic terms one obtains 



U{ft b , /i in ) = ( 1 + 2~ S rf) ^ ' \ +_ 4^ 2~ g ^ ' ' ^ ^ 



where 7 (0) is the coefficient of a s in the one loop anomalous dimension matrix that de- 
scribes the mixing of operators: 

? = i^ (0) + ■ ( 16 ) 

Note that it is ■y^ T and not that enters (1151) . Moreover, in the study of the /ij n 
dependence in the case of the scalar exchange one has to take into account that in this 
case the m 2 (/i in ) dependence is hidden in the coefficients (A % ^ R (H)) 2 . 

Considering then the cases of colourless gauge bosons and scalars we find that the following 
quantities should be /Xi n -independent: 

R Smge = U(fi b , fMin)C(^) , (17a) 

UT** = U(n b , /i in )C ? ( / u in )m 2 (/i in ) . (17b) 

For VLL (VRR) the column vector C(fi) is just a one-dimensional one, while it is two- 
dimensional for LR and SLL (SRR) systems. We write next 



C^ m )=Co-^p^K\og^, (18) 



where we suppressed /i in independent 0(a s ) terms. Moreover, the leading logarithm at 
0(a s ) in m 2 (fi in ) is given in 

mV in ) = m\M) (l + ^^ 0) log , (19) 

with 7w governing the scale dependence of the quark masses in QCD. 

Imposing (fTTj) . the conditions for K to ensure [i m independence of resulting amplitudes 
in these two cases read 

^gauge = JL^Co , (20a) 



K 



scalar 



>(0)T 

- + % 



7 ' ' -wi 



C , (20b) 



where 1 is a unit matrix. Thus the coefficients of logarithms in C(/j,- m ) can be found 
without the calculation of the loop diagrams in Fig. [T]l3]but the formulae in Eq. ( 120]) serve 
as a useful check of our results for logarithmic terms. These terms are renormalization 
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scheme independent and while cancelling the /ij n dependence of U /ii n ) in perturbation 
theory cannot remove its renormalization scheme dependence at the NLO level. To this 
end the 0(a s ) non-logarithmic terms have to be calculated which constitutes the main 
new result of our paper. 



In order to be able to use the relations in Eq. (120]) we recall the relevant one-loop anoma- 
lous dimension matrices [TU1: 



7 ( ° )VLL = 6 - | = 4 , 7 £> = 6C F = 3N - A = 8, (21) 



^=(f^%)Kt-W=6-»» (22) 



»_l -eiv + e + l j-i \_ /-2 7 ^ + e i-i \_ /-io | . 

i -24 -f 2iV + 6-|;-l -24-41 ^ 0) I + 6/"V-40#) (2d) 



Inserting these formulae into (|20p we indeed verify that the coefficients K gauge and K 8calaT 
of logarithmic terms calculated by us are correct: they are consistent with the /ij n - 
independence of the resulting physical amplitudes at this order of perturbation theory. 

It is instructive to compare the structure of the cancellation of the /z in -dependence in the 
case of the gauge boson exchange with the one of the scalar exchange: 

• In the gauge boson case the /independence of U (//&, /ii n ) can only be cancelled by 
the corrections calculated by us. 

• The case of scalar exchange with LR couplings is quite different. Here the //in- 
dependence of U(fib, A*in) is totally cancelled by the one of the m 2 (/i in ) so that 
even without our corrections the amplitudes are /i in independent. Indeed the C^f 
coefficients in the scalar case do not contain any logarithmic terms at 0(a s ). This 
type of cancellation can be traced to the fact that the anomalous dimension of 
the Q\ K operator is as seen in Eq. (I2"2l up to the sign equal twice the anomalous 
dimension of the mass operator. The role of our calculation is then the removal of 
the renormalization scheme dependence. 

• In the case of SLL operators the cancellation in question is not as pronounced 
because as seen in Eq. (123]) the anomalous dimensions of the relevant operators 
receive additional contributions beyond 7^- In this case our calculation removes 
both scale and renormalization scheme dependences. 
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5 Mixing Amplitudes at the NLO Level 



5.1 Preliminaries 



Having calculated 0(a s ) corrections to the Wilson coefficients at the matching scale /ii n 
we can obtain the complete NLO expressions for various tree-level contributions to the 
off-diagonal elements M 12 for the K and B s d systems. We will present only explicit 
expressions for AS* = 2 transition. Analogous expressions for B s _d systems can be easily 
obtained in the same manner. 

In presenting our results we will use the so-called P? QCD factors of pj] that include both 
hadronic matrix elements of contributing operators and renormalization group evolution 
from high energy to low energy scales. These factors depend on the system considered, 
depend on the high energy matching scale and depend on the renormalization scheme 
used to renormalize the operators. The formulae for these factors have been given in [TT] 
in the NDR-MS renormalization scheme of [22]. This scheme dependence is cancelled 
by the non-logarithmic 0(a s ) corrections calculated by us. The logarithmic corrections 
cancel the scale dependence of P? as explained in the previous section. 

The formulae for various contributions to M\2 are easily obtained from the Hamiltonians 
presented in Section |3] by replacing the operators by the corresponding hadronic matrix 
elements in Eq. (J2J) denoted there shortly by P^(fii n ). For completeness we recall the 
formulae for P^(fj, in ) [TT] : 

^i VLL (^in) = [vOml, A*i»)]vll S i VLL ( 24 ) 
P^ilMn) = ~\ [7/ n (/iL,^n)] LR [B^ R (fx L )] eff + ~ fo^, // in )] LR [5 2 l Vl)] cfr , (25) 

P 2 L V in ) = -i hi 2 (/x L ,/i in )] LR [B^ L )] eS + ~ M/^ta)] L R K R (/i L )] cfr , (26) 



A SLL (^i 



-r kl(/iL,/iin)] S LL [ 5 l LL (^)]efr - o [^l(/iL,/iin)] SLL [S 2 SLL (/1 L )] eff , (27) 



>SLL 



(^i 



"a [Wl(j*L, fMn)] SUj I 5 !^/^ 



eff 



2[^22(/iL,/iin)]sLL^2 LL (^)] eff , (28) 



where hl is a low energy scale at which hadronic matrix elements are evaluated. Explicit 
formulae for the QCD-NLO factors [^(/al, Hm)] a are given in [TT] . 

The effective parameters [-B"(/A£,)] e g are defined in the case of K° — K° mixing (fii = 
2GeV)by 



eff 



m K 



m a (ji L ) +m d (ii 1 



25.61 



98 MeV 



B!(»l). (29) 



In the case of B% s — B% s mixings one has to make the replacements jii —> and mj^F^ — > 
m 2 B Fj . Then in the case of B% — B® system (at fi b = 4.6) 



eff 



m B 



m b (fi b ) + m d (p b ) 



1.68 



4.08 GeV 



m b (fi b ) +m d (n b )_ 



B?W, (30) 
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with an analogous formula for the B s — B s system. 

We list now the final NLO expressions for the mixing amplitudes. 

5.2 Colourless gauge boson 



A«(A)A«(A) 



(31) 



The relevant Wilson coefficients are given in Section 13.2.11 
5.3 Colourless scalar 



M* 12 (AS = 2) = - (A ^ ))2 [C 1 SLL (^)P 1 SLL (^) + Cf L (/^)Pf L (M] 



K RR (^)pr R (^) + c 2 srr (^)p 2 srr (^)] (32) 



2M 2 H 



The relevant Wilson coefficients are given in Section 13.2.21 

We would like to emphasize that the formulae of this Section together with the QCD 
factors rjij presented in Section 3 of [llj and the coefficients Cj a (/ii n ) calculated in Section [3] 
of the present paper are valid for the tree level contributions of colourless bosons and 
scalars in any extension of the SM in which such contributions are present. The only 
model dependence enters through the couplings A l [ R and the gauge boson and scalar 
masses. In particular the coefficients Pf are universal in a given meson system and 
given renormalization scheme. With our normalization of Cf (/im) also these coefficients 
are universal except for the scale /ij n of NP and the same renormalization scheme used 
to evaluate Pf. In the process of multiplying Cj a (/i in ) and Pf terms 0(a 2 s ) have to be 
removed. 



6 Numerical Analysis 

We will now compute the size of 0(a s ) corrections and their impact on the reduction 
of the unphysical /x in -dependence present in the analyses in the literature. It should be 
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recalled that the actual size of the corrections calculated by us is not the most important 
result as these corrections are given in a particular renormalization scheme, the NDR-MS 
scheme. They could be different in another scheme. But then also the Pf factors would 
be different so that the final result for the physical mixing amplitudes M 12 would be 
renormalization scheme independent up to 0(a 2 ) corrections, that is NNLO corrections. 
Thus the important result of our paper is that we provide for the first time mixing 
amplitudes resulting from tree level decays including NLO QCD corrections that are 
renormalization scheme independent and which do not depend on a precise choice of the 
matching scale. Both statements are valid up to NNLO corrections. 

There are four linear combinations of P" and of the Cf that should be scale and renor- 
malization scheme independent. As we normalized the non-vanishing C? at LO to unity, 
these combinations are model independent. The model dependence enters only through 
the fermion-boson couplings and heavy boson masses characteristic for a given model. 
The linear combinations in question in the case of the K system are given as follows: 

Gauge Bosons: 



Rf = C 1 VLL (^)P 1 VLL (^a) , (33) 
Rf = C 1 lr (/. a )P 1 lr (^a) + C LR (MP 2 LR (M , (34) 

with the P" factors evaluated using hadronic K° — K° matrix elements at scale [it, — 
2 GeV. 

Scalars: 



P 



P 



x 



(cf^MA^GM + c 2 SLL (MPf LL (M) 
(c lr (/^)p lr (M + c* LR (/i„)p 2 LR (M) 



m 2 (M H ) 
m 2 (fj lH ) 



m 2 (M J 



H 



(35) 
(36) 



Note that even if the quantities P 2 and P4 appear at first sight to be the same they differ 
from each other because the Wilson coefficients of the operators for the gauge boson 
case are different than for the Higgs case. 

In Figs. Owe plot Rf as functions of the matching scales setting as an example the masses 
of gauge bosons and scalars to 1 TeV. Since the /i in dependence is the same in the B and 
K system we do not show the results for B mesons. They only differ in magnitude from 
each other because the hadronic matrix elements hidden in P" are different in these two 
meson sectors. Moreover, QCD effects are generally stronger in the K system because 
the renormalization group evolution is over a larger range of scales. 

In order to find numerical values of Pf one needs the values of the corresponding non- 
perturbative parameters Bf defined in [UJ . These are given in terms of the parameters 
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Figure 5: The quantities Rf defined in Eq. (Tff3j)-( Tff6j) as a function of (jLa,h for Ma,h = 
1 TeV. In the first row are the results for the gauge boson exchange and in the second row 
for scalar exchange. The dashed blue line represents the LO calculation for the Wilson 
coefficients, the dotted green line includes only logarithmic 0(a s ) corrections which reduces 
the scale dependence and the solid red line includes both logarithmic and non-logarithmic 
0(a s ) corrections that is important for the cancellation of scale and scheme dependence. 



Bi used in El [23J HQ as follows: 

^ LL (^ ) = B™(/io) = B 1 ( tM >) , (37a) 

B™-(jm>) = B b (jm>) , (37b) 

5 2 LR (/i ) = B 4 (//o) , (37c) 

Bf^M = 5 2 SRR (/io) = B 2 (fi ) , (37d) 

5 2 SLL (/i ) = S 2 SRR (/io) = \b 2 {^) - 2 -BM . (37e) 



The values for Bi in the MS-NDR scheme extracted from [23J for the K° — K° system are 
collected in Table HJ together with the relevant value of fj, . 

In each case we vary conservatively the matching scale between 300 GeV and 3000 GeV 
and show in each plot three curves: 
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B 1 


B 2 


B 3 


5 4 


#5 


A*0 


K°-K° 


0.57 


0.68 


1.10 


0.81 


0.56 


2.0 GeV 



Table 1: Values of the parameters Bi in the MS-NDR scheme obtained in [23] for K° — K° 
system. The scale /xq at which Bi are evaluated is given in the last column. 

• The result without the inclusion of 0(a s ) corrections as used until now in the 
literature (blue, dashed line). 

• The result after only logarithmic terms in the 0(a s ) have been included. They are 
crucial for the cancellation of the matching scale dependence (green, dotted line). 

• The result after the inclusion of non- logarithmic 0(a s ) corrections that are crucial 
for the cancellation of the renormalization scheme dependence (red, solid line). 

These plots are self-explanatory and we make only a few comments: 

• As expected the dotted green lines cross the dashed blue lines at /j,a,h = 1000 GeV. 

• The crossing point between the solid red and dashed blue lines is generally at [±a,h 
that differs from the mass Ma,h = 1000 GeV of the exchanged particle. The modest 
size of these differences shows that these corrections are small, of 0(5)% in the 
NDR-MS scheme. 

• In the case of gauge boson exchanges the matching scale dependence of roughly 
10(25)% for VLL(LR) operators in the range considered, when 0(a s ) corrections 
are not taken into account, has been reduced down to 1 — 2% after the inclusion of 
these corrections. 

• In the case of scalar exchanges the matching scale dependence in the case of the 
SLL sector of roughly 10% in the range considered, when 0(a s ) corrections are not 
taken into account, has been reduced down to 1 — 2% after the inclusion of 0(a s ) 
corrections. 

• In the case of LR operators in the scalar case there is basically no left-over scale 
dependence even in the LO as we explained at the end of Section HI 

This reduction of scale uncertainties cannot be appreciated at present in view of significant 
uncertainties in the values of the Bi parameters, as seen in Tabled], but the recent advances 
in lattice calculations allow for optimism and we expect that during this decade these 
uncertainties could be reduced below 5% and then the calculations presented here will 
turn out to be important. 
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7 Summary 

If there is a new physics at distance scales as short as 10~ 19 — 10~ 21 m, it will manifest 
itself primarily not through penguin and box diagrams as in the SM but through tree 
level FCNC processes. The best known examples of such NP are various versions of 
the so-called Z' models in which new neutral heavy weak boson (Z 1 ) mediate FCNC 
processes already at tree level. Gauged flavour models with new very heavy neutral gauge 
bosons and Left-Right symmetric models with heavy neutral scalars are other prominent 
examples where tree-level contributions to AF = 2 amplitudes are present. 

Effective tree-level contributions to AF = 2 observables can also be generated at one loop 
in models having GIM at the fundamental level and Minimal Flavour Violation of which 
Two-Higgs Doublet Models with and without supersymmetry are the best known exam- 
ples. In models with heavy vectorial fermions that mix with the standard chiral quarks 
and models in which Z° and SM neutral Higgs H° mix with new heavy gauge bosons and 
scalars in the process of electroweak symmetry breaking also tree-level contributions to 
AF = 2 processes mediated by Z° and SM neutral Higgs H° are possible. In all these 
extensions new local operators absent in the SM are generated having Wilson coefficients 
that are generally much stronger affected by renormalization group QCD effects than it 
is the case of the SM operators. 

Present studies of renormalization group QCD effects performed at the NLO level in many 
extensions of the SM use the so-called QCD factors [TTJ that include both hadronic 
matrix elements of contributing operators and renormalization group evolution from high 
energy to low energy scales. These factors represent the dominant part of any NLO 
QCD analysis but do not take into account 0(a s ) corrections to Wilson coefficients at 
the matching scale which separates the full and effective theories. Therefore basically 
all published calculations that considered tree level decays suffer from some unphysical 
scale and renormalization scheme dependences. While presently these unphysical effects 
are much smaller than the uncertainties in the hadronic matrix elements of contributing 
operators, the situation could change in this decade due to important progress in lattice 
simulations with dynamical fermions [251 [261 EZl [281 EHl [301 EU [32l [33]. 

While a general calculations of 0(a s ) corrections to Wilson coefficients, when the leading 
contributions come from loop diagrams, is very model dependent, a rather general analysis 
can be done for tree level exchanges so that the final results depend only on the couplings 
of exchanged bosons (vectors and scalars), on the QCD colour factors and the QCD 
coupling constant. 

The main goal of our analysis was to provide analytical formulae for arbitrary number 
of colours (N) for the 0(a s ) corrections in question in the case of tree level AF = 2 
processes mediated by heavy colourless gauge bosons and scalars. The results for the 
Wilson Coefficients and effective Hamiltonians for these cases are collected in Section El 
while the corresponding mixing amplitudes at the NLO level that combine our results with 
the known QCD factors are presented in Section [51 In Section [6] we demonstrated that 
the unphysical scale dependences have practically been removed by our calculations. This 
is particularly important for the LR system in the case of gauge boson exchanges, where 
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the scale dependence at LO is sizeable. The Appendices collect certain technicalities 
about the evanescent operators and intermediate results in the full and effective theories 
which should enable interested readers to check our calculations. 
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A The issue of Evanescent Operators 

It is well known that in the process of NLO calculations in the NDR-MS scheme, where 
ultraviolet divergences are regulated dimensionally, the so-called evanescent operators 
that vanish in D = 4 dimensions have to be considered. They arise in particular when 
complicated Dirac structures are projected onto the chosen basis of physical operators. 
The treatment of this operators in the process of matching considered by us must be 
consistent with the one used in the calculation of two-loop anomalous dimensions. 

We have used the QCD factors from [UJ which were based on the two-loop anomalous 
dimensions of operators calculated in [10]. Therefore it is mandatory for us to treat 
evanescent operators appearing in our calculations in the same manner as done in [TO] . 
Now, the latter paper used the treatment of evanescent operators as proposed in the 
context of the formulation of the NDR-MS scheme introduced in [22]. The virtue of this 
treatment is that the evanescent operators defined in this scheme influence only two-loop 
anomalous dimensions. By definition they do not contribute to the matching and to the 
finite corrections to the matrix elements of renormalized operators calculated by us. They 
are simply subtracted away in the process of renormalization. This issue is summarized 
in Section 6.9.4 of [20], where further references can be found. A very important paper 
in this context is [31]. Therefore effectively the calculations presented here were based on 
the projections listed in the next appendix that leave out the evanescent operators on the 
r.h.s. 

In this context we should warn the reader that the NDR-MS scheme used in [35], while 
sharing all the virtues of the scheme of [22] uses different projections in the SLL sector. 
This implies different two-loop anomalous dimensions for the SLL operators, that is also 
different P^ hh and also different 0(a s ) terms in Cf LL (/i in ) so that the physical amplitudes 
are the same in both schemes. The relation between these two schemes has been worked 
out in [12] . 

Another technical issue is related to the Fierz- vanishing evanescent operators, which have 
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to be considered when one wants to relate the operators with non-singlet structures like 

Qf L = (s a P L d^P L d a ), (38) 
Qf L = {s a a^P L d^{s^^P L d a ) (39) 

to the operators Qf L ^ used by us. In D = 4 — 2e dimensions the usual D = 4 identities 

or — ig? LL + ig2 LL , (40) 
Qf L D = 4 egf L + ^gf L . (41) 



do not work and one has to add evanescent operators on the r.h.s of these equations. 
However again, even if the inclusion of the latter operators was relevant for the calculation 
of two-loop anomalous dimensions of the physical operators in [10], it turns out that they 
do not contribute to the matching as long as the infrared divergences are not regulated 
dimensionally. As in our paper we regulated such divergences by a no n- vanishing p 2 , 
we can use the relations in Eq. ( 14*01 and Eq. (I41j) without taking the Fierz-vanishing 
evanescent operators in question into account. As discussed in [10] these "problems" are 
absent in the case of other operators. 

B Projections 

We list projections of all Dirac structures on physical operators that we encountered in our 
calculations. These projections correspond to the so-called "Greek method" as described 
in Section 6.9 of [22]. The evanescent operators, relevant in this renormalization scheme 
only at two-loop level, are defined as the differences of l.h.s and r.h.s of these equations. 
See [10] for more details. We define = | [y v , 



7a7/J7„(l ± 75hV ® 7 M (1 ± 75) = 4(1 - 2e) 7 „(1 ± 75) ® 7 M (1 ± 7s) (42a) 

7 M (1 ± 75)7,7^ ® y (1 ± 75)7 V = 4(4 - e) 7m (1 ± 75 ) <g) 7 "(1 ± 75) (42b) 

7/1 (l ± 75 ) 7q7 ^ ® 7 V7 M (1 ± 7s) = 4(1 - 2e) 7 „(1 ± 75) ® 7 M (1 ± 7s) (42c) 

7q7 ^(1 ± 75 )7V ® 7 M (1 T 75) = 4(1 - 2e) 7m (1 ± 75) ® 7 M (1 T 7s) (43a) 

7m (1 ± 75 )7 a7/3 <g> 7^(1 t 75)7 V = 4(1 + e) 7/i (l ± 7s) ® 7 M (1 T 7s) (43b) 

7 M (1 ± 7s)7a7/3 ® 7 /3 7 a 7"(l =F 7s) = 16(1 - e)j^(l ± 75) ® 7^(1 =F 75) (43c) 
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7^(1 T 75)7 V ® (1 ± 75) = 16(1 - e) (1 T 75) ® (1 ± 75) (44a) 
(1 =F 75)7m7, ® (1 ± 75)7^ = 4(1 + e) (1 =F 7s ) <g> (1 ± 75) (44b) 
(1 T 75)7m7, ® 7^(1 ± 75) = 4(1 - 2e) (1 =F 75) ® (1 ± 7s) (44c) 



7,7^(1 ± 75)7^ ® (1 ± 75) = 16(1 - e) (1 ± 75) ® (1 ± 75) 

(1 ± 75)7^ ® (1 ± 75)7^7" 

= (4 - 2e) (1 ± 75) ® (1 ± 75) + <v(l ± 75) ® ^(l ± 75) 

(1 ± 7s)7m7. ® 7V(1 ± 75) 

= (4 - 2e) (1 ± 75) ® (1 ± 75) - <V(1 ± 75) ® ^(1 ± 75) 



(45a) 
(45b) 

(45c) 



rfrf a V»Q ± 75 ) 7/37q ® ± 75) = 

<V(1 ± 75 )7a7/3 ® ^(1 ± 7 5 hV 

= (48 - 80e) (1 ± 75 ) ® (1 ± 75 ) + (12 - 6e) <v(l ± 75) ® <7^(1 ± 75) 

<V(1 ± 75 )7a7/3 ® 7 /3 7 a ^(l ± 75) 

= -(48 - 80e) (1 ± 75 ) ® (1 ± 75 ) + (12 - 14e) <v(l ± 75) ® ^(1 ± 75) 

C Matrix Elements of Operators 

After quark wave function renormalization and operator renormalization we get 

■ 2 r ~ 



(Q\ R ) = 
(Q2 R ) = 



1 , on a * 1 ^ 1 q s A A 1 ' 



1 + 2C F ^- log 

47T — £T 



1 + 8C F -±- ( 1 + log 

47T 



— p 2 



a, / 1 



q7 ll 

2N 1 2 



iLR 



1 



LR 



(46a) 
(46b) 

(46c) 



(47) 
(48) 
(49) 
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(Qi LL ) 



a 



1 + SC F -^- 1 + log 



47T 



Arc 



2 + log 



— 



— p 2 
N 



Qi 



SLL 



Qf L + 3Q* 



SLL 



48^flo~ " 

47T 



4iV 
2 + 3 



p 



AN 



N -2 



Qf L + 3Q? 



4tt V AN 



SLL 



(50) 



iV + 2 Qf L + ^Qf L 



(51) 



We remark that for the matching performed in this paper only the 0(a s ) corrections to 
the matrix elements (Q7 LL ), (Qi R ) and (Qf LL ) matter. We give the remaining matrix 
elements as they are relevant for NLO matching when coloured gauge bosons and scalars 
are exchanged [19J. 



D Amplitudes in the Full Theory 

Colourless gauge boson exchange (after quark wave function renormalization): 



(A\ 



VLL\gaug- 



,_ (Af(A)) 
2M\ 



2C, 



Air 



log 



/i 2 N- 1 a s (\ 



-p z 



N 4tt \2 



+ log 



Ml 
—p 2 



Qi (52) 



^LR^ gauge _ A S l(A)Ax(A) 



M 2 A 



An 



~P 



a 



1 , M 2 A 
- + log^f 
An \ 2 —p z 



3 

N 



(53) 



QLR + ggLR 



Colourless scalar exchange (after quark wave function and quark mass renormalization): 



^^LR^scalar 
^^SLL ^ scalar 



Af(H)A^(H) 



M 2 H 



(Agf(ff))' 
2M 2 H 



1 + 8C ^( 1 + l0 ^ 



a 



SC^fl + lo, 



a s fl , M 2 
A-,{2 +h ^ 



P 

N-2 



AN 



Qf LL 

Ql LL + 3Qf L 



(54) 



(55) 
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